o 
o 
o 

(N 

m 






Transport Properties and Density of States of Quantum Wires with 
Off-diagonal Disorder 

P. W. Brouwer ab * , Christopher Mudry act and Akira Furusaki d * 

a Lyman Laboratory of Physics, Harvard University Cambridge, MA 02138, USA 

b Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY 14853, USA 

c Paul Scherrer Institut, CH-5232, Villigen PSI, Switzerland 

d Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 



Ctf 



O 
o 



Abstract 

We review recent work on the random hopping problem in a quasi-one-dimensional geometry of N coupled 
chains (quantum wire with off-diagonal disorder). Both density of states and conductance show a remarkable 
dependence on the parity of N. The theory is compared to numerical simulations. 
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1. Introduction and Results 

With the realization that the problem of An- 
derson localization is amenable to a RG analy- 
sis came the understanding that some transport 
and spectral properties of a quantum particle sub- 
jected to a weak random potential depend qual- 
itatively only on dimensionality and the symme- 
tries of the random potential [jjj . For essentially 
all types of randomness the dimensionality two 
of space plays the role of the lower critical di- 
mension: Below two dimensions and for arbitrary 
weak disorder a quantum particle is always lo- 
calized (i.e., the wavefunction is insensitive to a 
change in boundary conditions, or, equivalently, 
exponential decay of the conductance with sys- 
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tern size), whereas a finite amount of disorder is 
needed to localize a quantum particle in three di- 
mensions. 

A notorious exception to this rule is the case of 
a particle on a single chain with random nearest 
neighbor hopping which appears in many rein- 
carnations, e.g., random XY spin chains p| and 
diffusion in random environments S . In all these 
cases, the cause underlying the different behavior 
of random systems with off-diagonal disorder is 
the existence of a sublattice, or chiral, symmetry, 
which is absent for diagonal disorder |jj . 

The purpose of this contribution is to review 
our recent work (together with Simons and Al- 
tland fa]) on the random hopping problem in a 
quasi-one-dimensional "wire" geometry (&-p| . We 
consider the density of states (DoS) and con- 
ductance, and focus on the differences between 
the cases of off-diagonal and diagonal disorder 
in the localized regime. (For off-diagonal dis- 
order, we assume absence of diagonal disorder.) 
The quantum wire is the logical intermediate be- 



tween one and two dimensions. [J The quasi-one- 
dimensional geometry allows us to find exact so- 
lutions for the conductance at the band center 
e = and the DoS near e — 0, while, in con- 
trast to purely one-dimensional counterparts, it 
still shows a crossover from a diffusive to a local- 
ized regime, as is the case in two dimensions. 

In the diffusive regime, differences between off- 
diagonal and diagonal disorder are limited to 
small quantum corrections to the DoS and con- 
ductance. The DoS near e = follows from the 
so-called "chiral" random matrix theory (RMT), 
which was originally proposed in the context of 
quenched approximations to the QCD Hamilto- 
nian O). In chiral RMT, all eigenvalues come in 
pairs ±e. For N coupled random hopping chains, 
level repulsion between the smallest eigenvalue 
above the band center e = and its mirror image 
then leads to a suppression of the DoS near e = 0, 



dAf/de (xv^ie/A) 13 - 1 , < e < A. 



(1) 



Here dAf/de is the DoS per unit volume, vf is the 
Fermi velocity, A = 2-kvf/NL is the mean level 
spacing, and the symmetry index f3 = 1 in the 
presence of both time-reversal and spin- rotational 
symmetry and j3 = 2 (4) if time-reversal symme- 
try (spin-rotational symmetry) is broken. Several 
level spacings away from the center of the band, 
the mirror symmetry of the spectrum plays no im- 
portant role, and the level statistics were found to 
agree with those of the standard RMT. Similarly, 
different quantum corrections to the conductance 
are obtained for energies close to the band center 

On the other hand, beyond the diffusive regime, 
i.e., for sample lengths L > £, where £ is a 
crossover length scale (to be defined after Eq. 
(0) below), quantum corrections are dominant, 
and the differences between diagonal and off- 
diagonal disorder can be much more dramatic. 
While for diagonal disorder, the DoS is constant 
as e — ► 0, Dyson showed that for off-diagonal dis- 
order dAf/de diverges upon approaching the band 



1 Two dimensional models with off-diagonal disorder are 
also studied in the context of strongly interacting electron 
systems |9|. We refer the reader to Ref. [ho(| for a review 
of work done on the two dimensional case. 



center JL2| : 

dAf 
de ti~i- 3 



< eijvF < 1. 



(2) 
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Note that, in contrast to the diffusive regime (||), 
the energy scale that governs the divergence does 
not depend on the system size, nor does the form 
of the singularity depend on the symmetry index 
(3. For N coupled chains, Dyson's result (||), re- 
mains valid for odd Nq but not for even N ||: 
For an even number of chains, the DoS diverges 
only logarithmically for (3 = 1, whereas dAf/de 
shows a pseudogap for (3 = 2, 4, 



dAf 1 
de vf 



0-1 



V F 



0< — <1. (3) 

V F 



The relevant energy scale, however, remains the 
same as in Eq. (^). The even-odd effect in the 
DoS around e = is accompanied by an even-odd 
effect for the conductance g at the band center [&- 
M|. For odd N, there is no exponential localiza- 
tion; g has a broad distribution, with fluctuations 
that are of the same order as the average, 



(In g) 



a L , 8 7T-2X 

- 4 V7TT' varlng= ' , 4 



where the crossover length scale £ = (f3N + 2 — 
0)1/ j3, and £ is the mean free path. (For small N 
there are corrections due to the appearance of a 
second dimensionless parameter needed to char- 
acterize the off-diagonal disorder, see Ref. M for 
details.) In contrast, for even N, g is exponen- 
tially small, its distribution being close to log- 
normal, 



(In g) 




%L 8(tt-1)L 
— -, var In g = — — .(5) 



In this case, £ can be interpreted as the local- 
ization length. Away from the band center and 
for diagonal disorder, (lng) = — (l/2)var lng = 
-2L/£, with £' = {f3N + 2 - [3)1, irrespective of 
the parity of N , as for quantum wires with diag- 
onal disorder jR|. Note that for even N, despite 
the similarities, differences between off-diagonal 

2 The proportionality constants in Eqs. (H) and (fa) may 
depend on TV. 



and diagonal disorder persist, though they are 
more subtle. 

In the remaining sections of this paper, we re- 
view the theory describing these two parity ef- 
fects. 

2. Two microscopic models 

2.1. Lattice model 

The random hopping model is a lattice model 
with nearest neighbor hopping only, where the 
hopping amplitude takes a different (random) 
value between adjacent sites. In our case we con- 
sider a two-dimensional square lattice, for which 
the Schroedinger Equation takes the form 



W>i. 



2 , {ti,j:i±l,j' l Pi±l,j + ti,j;i,j±l1pi,j±l) ■ (6) 



Hermiticity requires ti t j-i±i t j = (ti±i,j;i,j)* an d 
tij-ij±i = {ti.j±i;ij)*. We consider a quasi-one- 
dimensional geometry (length L of the disordered 
region much larger than its width Na, where N is 
the number of chains and a the lattice constant), 
with open boundary conditions in the transverse 
direction: to,j;i,j = tN,j;N+i,j = 0. On the left 
(j < 0) and right (j > L/a), the disordered re- 
gion is attached to ideal leads (tjjjiiij = ill, 
tij. t ij±i = t±. for j < or j > L/a). In the 
disordered region, the hopping amplitudes show 
fluctuations around the average values i|| and tj_, 
respectively. The fluctuations are real, complex, 
or quaternion, for the symmetry classes (3 = 1, 2, 
and 4, respectively. 

The random flux model is a special version of 
the random hopping model for which each plaque- 
tte is threaded by a random flux. In that case, 
the hopping amplitudes have magnitude t and a 
random phase 4>ij-i±i t j, so that the total phase 
accumulated going around a plaquette equals the 
phase through the plaquette in units of the flux 
quantum $o. 

The random hopping model is special because 
of the existence of a sublattice symmetry: Under 
a transformation 

iHj - (-l) i+i ^,i, 

the Hamiltonian changes sign. Hence, if e is an 
eigenvalue of H, then — e is so as well. The band 



center e = is a special energy with respect to 
this transformation. As we shall see below, it is 
the existence of this extra symmetry that causes 
the spectral and transport properties of the ran- 
dom hopping model at the band center to be so 
dramatically different from those of models with 
on-site disorder. For historical reasons, the sub- 
lattice symmetry is referred to as chiral symme- 
try. 

2.2. Continuum model 

While the lattice version (||) of the random hop- 
ping model is the version that is usually consid- 
ered in the literature, we use a different (contin- 
uum) model for our analytical calculations, 

H cont ip(y) = [ia 3 d y + cr 3 v(y) + a 2 w(y)} if){y). (7) 

Here rp is a 2N component vector (elements of ij; 
occur in pairs that correspond to left and right 
movers), v and w are N x N random Hcrmitean 
matrices, and the a^ (/i = 1,2,3) are the Pauli 
matrices. In Eq. fm) and below we choose our 
units such that the Fermi velocity vf = 1. The 
chiral symmetry is now represented by 

Cl7i con t<Ti = — Ticont- 

For (3 = 1, w (v) is (anti-) symmetric; for (3 = 4, 
w (v) is (anti)-self dual. The disorder potentials 
v and w are chosen independently and Gaussian 
distributed with zero mean and with variance 
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(w lj (y)wl l (y')} = 



S(y - y') 

5(y - y') 



\ s ik s ji ~ ^ 5 u s jk 



($ik s ji 
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respectively. (For small N a slightly more gen- 
eral form of the variance is needed, see 0.) The 
justification for this change of models is our an- 
ticipation that, for weak disorder and for a suf- 
ficiently large system size, only the fundamental 
symmetries of the microscopic model are relevant. 
In order to verify this assumption of universality, 
we compare the theoretical predictions based on 
Eq. (0) below to numerical simulations of the lat- 
tice model (g). 



3. Fokker-Planck approach 

3.1. General method 

Spectral and transport properties of the con- 
tinuum model (m) can be studied in a unified way 
from the statistical distribution of the reflection 
matrix r(e) of the disordered wire. The reflection 
matrix is defined from the relation between the 
amplitudes ip 1 ^ and ijj° h of an incoming and an 
outgoing wavefunction at energy e on the same 
(left) side of the disordered region (Fig. 1), 



v>; 



oL 



"(£)€ 



(8) 



Transport properties (the conductance G) follow 
from the eigenvalues of r^r, when the right side 
of the quantum wire is attached to a lead (Fig. 
la), 



r 2e2 



tr 1 



r'r) 



(9) 



while the DoS is obtained from the eigenphases 
of r when the quantum wire is closed at the right 
side |1 (Fig. lb) 



dN 1 / + d 

de 2ttiNL \ de 



(10) 



The distribution of r is computed using a RG 
approach: One calculates how the eigenvalues of 
r^r (in the case of the conductance) or the eigen- 
phases of r (in the case of the DoS) change when 
a thin slice is added to the disordered region (Fig. 
g) . If we place the thin slice to the left of the dis- 
ordered region, the change of the reflection matrix 
is given by 



-i+*i(l-rri) \t 1) 



(11) 



where r\ and r[ (t\ and t'i) are the reflection 
(transmission) matrices of the thin slice (Fig. ||). 
If the length SL of the thin slice is much smaller 
than a mean free path, these can be computed 
using the Born approximation, 

r x = -W +i(VW -WV), 

t x = l + iV-\{V 2 +W 2 )+ie5L, 



W 



-{VW -WV), 



t[ = l-iV-~{V 2 +W 2 )+i£5L, 



(12) 



where V = J dyv(y), and W — L dyw(y). 
Here we neglected terms that are of order 5L 2 . 
Together with the distribution of the disorder po- 
tentials w and v, Eqs. (O)-(O) define how the 
distribution of the reflection matrix r evolves with 
the length L of the quantum wire. The RG ap- 
proach can be represented in terms of a Fokker- 
Planck equation describing the "Brownian mo- 
tion" of the eigenvalues of r*r or eigenphases of r 
in the geometry of Fig , la and lb, respectively. 

Equations (jll| )-(|l2|) are exact for the contin- 
uum model (0)7 A different choice for the dis- 
tribution of v and w, or use of the lattice model 
@ instead of Eq. (j7|), would have led to differ- 
ent statistical properties of the scattering matrix 
for a thin slice. However, as we have verified nu- 
merically, such differences are irrelevant in the 
RG sense, i.e., they disappear for sufficiently long 
wires (longer than the mean free path £) and weak 
disorder (I much larger than the lattice constant 
a). 

3.2. DoS 

Since the quantum wire is closed at one end, 
the initial condition at L = for the RG equa- 
tion ([Fy) is r = 1. The RG flow ensures that 
the reflection matrix r remains a N x N unitary 
matrix for all L. Hence, for all L there exists a 
Hcrmitean N x N matrix $ with r = exp(i<I>). 
The eigenvalues <j>j of $ are the eigenphases of r. 
The change of $ under the increment ([ll|)-(|l2]) 
is 



cotl- 



,-(W+iV) , 



cot [ h eSL I e 



,-OV-iV) 



up to second order in V and W and first order in 
e. Taking the length L as a fictitious "time", the 
eigenphases (f>j of $ perform a Brownian motion 
on the unit circle, which is such that upon increas- 
ing L, they move (on average) counterclockwise. 
Integration of Eq. (|l(j) yields a relation between 
the <j)j and the DoS valid in the limit of large L, 



Af(s') = 



de 



dN 
de 



1 



a d(j)j 
2ttN ^ ~dL 



(13) 



[Note that for e — 0, <p ■ = for all L.\ Hence, to 
find the (average) DoS per unit volume, we com- 
pute the (average) "current' of the eigenphases 



cj)j moving around the unit circle and differenti- 
ate to energy. We remark that, in the absence of 
disorder, the angles 4>j move around at a constant 
speed oc e, resulting in a constant DoS. With dis- 
order, their motion acquires a random (Brown- 
ian) component, which dramatically affects their 
average speed, and hence the DoS. 

For a quantitative description a different pa- 
rameterization of the eigenphases of the reflection 
matrix is more convenient, 



tan(0j/2) = cxp(uj), 



(14) 



where the Uj are restricted to the two branches 
Imuj — and Imuj = ir in the complex plane 
(Fig. y). Starting from the RG equation (pj), 
one then obtains that the joint probability distri- 
bution P(ui, . . . , un) for the Uj obeys 



dP 



3 



d 



4 d 
-zzJ-^—J -2ecosh(u i ) 



J =Y[smh fj [(u j -u k )/2]. 



P 



(15) 



3<k 



This Fokker-Planck equation describes the mo- 
tion of N fictitious Brownian "particles" with 
coordinates Uj (j — 1, . . . , N) and diffusion co- 
efficient 4//3£, subject to a driving force F e = 
2e coshuj and a long-range repulsive two-body in- 
teraction F int — (2/£) coth[(uj — u k )/2}. The pa- 
rameterization (Um is such that a Brownian par- 
ticle with coordinate Uj that vanishes on one of 
the branches at ±oo reappears at the opposite 
branch (dotted arrows in Fig. 0) . 

Analysis of the Brownian motion described by 
Eq. ( p"5| ) then leads to the asymptotes (0J3) for 
the DoS in the random hopping model. A de- 
tailed derivation of these results based on an es- 
timate of the steady-state-current supported by 
Eq. ( |l5| ) can be found in Ref. || . Here we present 
only a brief sketch, focusing on the origin of the 
even-odd effect: The competition between the 
driving force F e and diffusion on the one side, 
and the long-range repulsive two-body interac- 
tion F; nt on the other side. The driving force 
F e pushes the particles to the right (left) on the 
lower (upper) branch and thus causes the nonzero 
steady-state-current. Diffusion keeps the parti- 
cles mobile where the driving force F s is small, 



thus enhancing the current, and hence the DoS 
||. The long-range repulsion F int , on the other 
hand, traps an equal number of the particles in 
two "traps" near u = ln(e£) and u = — ln(e£)+7ri, 
and thus tends to suppress the DoS. For even N, 
all particles are trapped {N/2 particles in each 
trap). Rare events caused by particles that are 
"thermally" excited out of their traps still allows 
for a small "current" , resulting in the small DoS 
in Eq. (||) (Fig. ||a). For odd N, however, only 
{N — l)/2 particles are trapped near each end 
point, while one particle is left alone, unaffected 
by the long-range repulsion (Fig. Hb). Motion of 
this particle is dominated by diffusion, which ex- 
plains the enhancement of the DoS for odd N. 

The agreement between the theory and nu- 
merical simulation of the lattice random hopping 
model is excellent as is depicted in Fig. 0. 

3.3. Conductance 

The chiral symmetry results in the symmetry 
relation r(e) = r{— e)^ M. Hence, at the band 
center, the reflection matrix is Hermitean. We pa- 
rameterize the eigenvalues of r as tanh Xj , where 
the Xj are N real numbers, and formulate the RG 
approach of Sec. 3.1 as a Brownian motion pro- 
cess for the Xj. Away from the band center, only 
the eigenvalues tanh Xj of the product r*r can 
be studied in the RG approach: In this case, the 
sign of the Xj has no relevance, since r'r does not 
change under Xj — > — Xj. 

In terms of the x 



3' 



the conductance reads 



= /J(l — tanh 2 Xj) — 2_] cosh 



(16) 



The RG equation ( |TT| ) then leads to a Fokker- 
Planck equation for the L-evolution of the joint 
probability distribution P{xi, . . . ,Xn) at the 
band center 



dP_ _ 1 y^ 8 

3 = 1 J 

J = I I | sinh(xj - 



x k )f. 



(17) 



Sufficiently far away from the band center, the L- 
evolution of the Xj 's is the same as for a quantum 



wire with on-site disorder, which is given by the 
well-known DMPK equation Jlgl 



dP__ 1 J^ 8 



dL 2£ 



j=i 



dxj 



J ^ J - 1P ) 



(18) 



J = Y[ | sinh 2 Xj - sinh 2 ^fc|^ J^[ | sinh(2:Ej)|, 

k>j k 

where £' = (/3JV + 2 - /?)£. 

Both Fokker-Planck equations describe the 
Brownian motion of iV particles with the coordi- 
nate Xj on the real axis and interacting through 
a hardcore repulsive two-body potential (Fig. a). 
Far away from the band center, a particle at po- 
sition Xj is also repelled by the mirror image at 
position — 
r^r under 



as is required by the symmetry of 
— > —Xj. For each Xj, the repulsion 



from the mirror image — Xj ensures that all Xj 
increase linearly with L, thus causing the expo- 
nential decay of g with L. On the other hand, at 
the band center the sign of Xj matters since there 
is no repulsion from mirror images. For odd N 
there will be no net force on the particle with co- 



ordinate X(j\^ 



L )/ 2 (for large L). Hence x 



(N+l)/2 



remains close to the origin, while all other Xj in- 
crease (or decrease) linearly with L. The coor- 
dinate X(jv+i)/2 dominates the conductance, its 
fluctuations leading to the (L/^) 1 ' 2 dependence 
of In g in Eq. (0). For even N no such "excep- 
tional" particle exists; all Xj increase (if they are 
positive) or decrease (if the are negative) linearly 
with L, causing the conductance to be exponen- 
tially small, cf. Eq. (||). 

We refer the reader to Refs. || for the detailed 
calculation of the moments of the conductance. 
The agreement between the theory and numerical 
simulations of the random flux case is excellent as 
is illustrated in Fig. ||. 
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FIG. pj. Boundary conditions used for the compu- 
tation of the conductance (a) and DoS (b). 
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FIG. El A thin disordered slice of length SL is 

added to the left of a disordered wire of length L. 
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FIG. W. Density of states, from numerical simula- 
tions for N = 1, 3 (a) and for AT = 2,4 (b). The data 
shown on the linear scale are computed for L — 200 
(L = 500 for N — 1) while the data in the insets are 
for L = 10 5 (L = 10 B for N = 1). For /3 = 1 and also 
for N — 1 the hopping amplitudes are taken from a 
uniform distribution in the interval [0.5,1.5], while for 
j3 = 2 the random flux model H is used, where the 
randomness is introduced only via the random phases 
of the hopping amplitudes. Results of an average over 
4 x 10 4 -10 6 disorder realizations are shown. 
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FIG. |3|. The L-dependence of the eigenphases 

of r is described in terms of a Brownian motion of 
fictitious particles with coordinate Uj where Im u = 
or n. If N is even, a repulsive interaction (see text) 
traps all particles near u — — In e£ or In e£ + i7r and 
prohibits them from moving around (a) . If iV is odd, 
one particle can diffuse freely around the branches, 
and thus leads to a significantly higher DoS than for 
even N (b). 
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FIG. H. For a wire with off-diagonal disorder, 

the L-dependence of the eigenvalues tanh 2 Xj of r'r is 
described in terms of a Brownian motion of fictitious 
particles with coordinate Xj on the real axis. If iV is 
even, all Xj (full circles) repel away from (a), while 
i(jv+i)/2 remains close to for odd N (b). In the 
case of diagonal disorder, all Xj repel from due to 
repulsion from mirror images (open circles) (c). 
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FIG. |6[ Average and variance of In g versus L/£ for 
N = 1,3 (a) and AT = 2, 4 (b). The hopping am- 
plitude t|| is taken from a uniform distribution in the 
interval [0.9,1.1], while t± is a real random numbers in 
[-0.2,0.2] for (3 = 1 and a complex random numbers 
with modulus < 0.1 for f3 = 2. The characteristic 
length ^ is obtained from a comparison with Eqs. (ft) 
and (pf). Results of an average over more than 2 x 10 4 
disorder realizations are shown. 



